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coefficients of the azimuthal distribution of the produced particles v\ to V3 within hybrid transport 
plus hydrodynamics approaches are summarized. 
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1. Introduction 


In heavy ion collisions at ultra-relativistic energies nuclear matter can be studied under ex¬ 
treme conditions. At zero baryo-chemical potential and finite temperature, a cross-over between 
the hadron gas at low temperatures and the quark gluon plasma at high temperatures has been 
established by lattice QCD [1,2] and also been confirmed by experimental measurements at the 
Large Hadron Collider (LHC) and the Relativistic Heavy Ion Collider (RHIC) in combination with 
dynamical models [3]. At finite net baryon densities, current lattice calculations are not applicable 
and insights about the structure of the QCD phase diagram can only be achieved by experimental 
exploration in combination with detailed dynamical modeling and phenomenological approaches 
providing input on the equation of state in the whole temperature (r)-baryo-chemical potential (Hb) 
plane. 

The previous and currently running programs at CERN-SPS and the recent beam energy scan 
(BES) program at RHIC provide a comprehensive data set on general bulk observables and selected 
fluctuation/correlation measurements. In the future, the second generation of experiments with im¬ 
proved capabilities for rare probes will meticulously explore the exciting energy regime accessible 
at EAIR, NICA, and stage 2 of the BES@RHIC. The second goal, besides obtaining insight on the 
type of the phase transition of strongly interacting matter, is to understand the properties of the 
state of matter that is formed, the QGP. 

One of the most promising observables is the anisotropic flow, that is measured in terms of 
Eourier coefficients of the azimuthal distribution of produced particles. Directed flow vi and elliptic 
flow V 2 are analysed extensively to gain insights on the transport properties and the geometry of 
heavy ion collisions. Within the last 5 years the whole plethora of odd higher coefficients has been 
studied to investigate the initial state and its fluctuations in more detail. The flow coefficients - as 
indicated by the name - are a sign of collective behaviour and quantify the response of the system 
to spatial anisotropies, translated to momentum space by pressure gradients. 

To connect the final state particle distributions to the quantities of interest, e.g. the equation 
of state or the viscosity to entropy ratio of QCD matter, detailed dynamical models are necessary. 
The current state of the art is to use a combination of 3-1-1 dimensional (viscous) hydrodynamic 
evolution and hadronic transport for the non-equilibrium evolution in the late stages of the reaction. 
There are various ways to model/parametrize the initial state for these calculation ranging from 
Glauber/CGC type models, to classical Yang-Mills evolution, AdS/CET or transport approaches. 
At lower beam energies the initial non-equilibrium evolution takes a non-negligible time and one 
needs to certainly pay attention to the type of initial dynamics that is required. 

In the present study, the calculations are performed within the Ultra-relativistic Quantum 
Molecular Dynamics (UrQMD) approach including a 3-1-1 dimensional ideal hydrodynamic evo¬ 
lution where appropriate. The main ingredients of the model are described in Section 2. After that. 
Section 3 contains a summary of the expectations for the beam energy dependence of anisotropic 
flow in heavy ion collisions and the thr'ee subsections 3.1, 3.2 and 3.3 are devoted to vi to V 3 in 
more detail. The last Section 4 summarizes the main conclusions. 
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2. Hybrid Approach at non zero jUg 

Heavy ion collisions are processes that do not create a system at fixed temperature and density 
that can be studied over an extended period of time, but rather span a fluctuating finite region in 
the phase diagram evolving dynamically. Hybrid approaches, based on (viscous) hydrodynamics 
for the hot and dense stage of the evolution, and non-equilibrium transport for the dilute stages of 
the reaction are successfully applied to describe bulk observables at high RHIC and LHC energies. 
These approaches combine the advantages of hydrodynamics and transport and apply both approx¬ 
imations within their respective regions of validity. The equation of state and transport coefficients 
of hot and dense QCD matter are direct inputs to the hydrodynamic equations, which allows for a 
controlled modeling of the phase transition. Microscopic transport on the other hand describes all 
particles and their interactions and provides the full phase-space distribution of produced particles 
in the final state, which allows for an apples-to-apples comparison to experimental data. 

At higher beam energies the standard picture of the dynamical evolution of a heavy ion re¬ 
action consists of non-equilibrium initial state dynamics providing the initial state for the viscous 
hydrodynamics equations and a hadronic afterburner to account for the separation of chemical and 
kinetic freeze-out. At lower beam energies, it is not clear how close the produced system really gets 
to local equilibration and if it thermalizes in the whole phase-space. In general, it is expected that 
dissipative effects are larger for low beam energies and the hadronic interactions gain importance. 
In addition, the widely used assumption of boost invariance breaks down and a 3-1-1 dimensional 
calculation is inevitable. Also, the equation of state and the transport coefficients are required not 
only as a function of temperature, but in the whole T — jis - plane. The finite net baryon current 
needs to be conserved during the evolution. Going beyond the standard hybrid approach by ex¬ 
ploring the non-equilibrium dynamics in the vicinity of a critical endpoint or a first order phase 
transition is the ultimate goal, where exploratory work has been performed so far [4]. 

The dynamical approach that is used to calculate the results presented below is the Ultra- 
relativistic Quantum Molecular Dynamics (UrQMD) approach coupled to ideal relativistic fluid 
dynamics for the hot and dense stage (SHASTA) [5, 6, 7]. This approach has recently been im¬ 
proved to include the finite shear viscosity during the hydrodynamic evolution as presented in an¬ 
other contribution to this conference and in [8]. The initial state is produced by generating nucleons 
according to Woods-Saxon distributions and computing dynamically the first binary interactions 
until the two nuclei have geometrically passed through each other. Event -by-event fluctuations of 
the positions of the nucleons and of the energy deposition per collision are naturally included. At 
4tart = ^ the energy, momentum and net baryon density distributions are calculated by represent¬ 
ing the individual particles with a three dimensional Gaussian distribution. Based on these initial 
conditions including fluctuating velocity profiles the ideal relativistic hydrodynamic equations are 
solved. There are different options for the equation of state, e.g. a bag model equation of state 
matched to a hadron gas at low temperatures to explore the sensitivity of observables to a strong 
first order phase transition (BM) or a more realistic equation of state provided by a chiral model 
that is fitted to lattice QCD at zero jJ-s, reproduces the nuclear ground state properties and incor¬ 
porates constraints from neutron star properties. At a constant energy density a hypersurface is 
constructed and particles are sampled according to the Cooper-Frye formula. The hadronic rescat¬ 
tering and decays are treated by the hadronic transport approach. This approach incorporates most 
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of the above mentioned ingredients that are necessary for a realistic dynamical description of heavy 
ion reactions at lower beam energies. The current strategy can be seen as a first step to take the 
well-established picture at high energies and explore how well it works at lower beam energies. 


3. Anisotropic Flow Observables 


charged particles, lykO.l 




Figure 1: (Color online) Left: The calculated beam energy excitation function of elliptic flow of charged 
particles in AuH-Au/Pb-rPb collisions in mid-central collisions (b=5-9 fm) with |y| < 0. l(full line). This curve 
is compared to data from different experiments for mid-central collisions (see [12] for refs). The dotted line 
in the low energy regime depicts UrQMD calculations with the mean field [9]. Fig. from [10]. Right: V 2 /e 
as a function of {l/S)dNch/dy for different energies and centralities for Pb-tPh/Au-tAu collisions compared 
to data [14]. The results from mid-central collisions (b=5-9 fm) calculated within the hybrid model with 
different freeze-out transitions and different definitions of the eccentricity are shown by black lines. The 
green full lines correspond to the previously calculated hydrodynamic limits [11]. Fig. taken from [12]. 

Anisolropic flow is the colleclive response lo geometrical structures in the initial state distribu¬ 
tion of heavy ion collisions. If there are strong enough interactions spatial anisotropies are getting 
translated to momentum space throughout the evolution. Therefore, anisotropic flow is supposed 
to be very sensitive to the equation of state and to the transport properties of the produced matter. 
It is quantified by the Fourier coefficients of the azimuthal distribution of the produced particles 

v„ = cos(n0-'T„) (3.1) 

where is the event plane and n specifies the coefficient of interest. There are many different 
two-particle and many-particle methods to measure the anisotropic flow coefficients with different 
sensitivities to non-flow and fluctuations. In this article, either the event plane method (some V 2 
results and V 3 ) or the theoretical definition with respect to the known reaction plane (vi = ( 2 ) 

P.x ' Py 

2 __ 2 

and V 2 = ( ^21 ^ 2 )) have been applied to calculate the flow observables. 

Px~^ Py 
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Figure 2: (Color online) Left: Fraction of elliptic flow generated during hydrodynamic evolution as a func¬ 
tion of beam energy. Right: Magnitude of V 2 {EP} at the beginning of hydrodynamical evolution (squares), 
immediately after particlization (diamonds) and after the full simulation (circles) in midcentral collisions. 
Fig. from [13]. 

Elliptic flow is one of the main observables to support the claim that the quark gluon plasma 
behaves like a nearly perfect liquid. Therefore, the beam energy dependence of elliptic flow is sup¬ 
posed to be sensitive to the phase transition and the changing transport properties at lower collision 
energies. Triangular flow is the first odd flow component that is mainly sensitive to fluctuations and 
would be zero, if the initial state is taken as average over many events. Therefore, the hope is that 
measuring the beam energy dependence of triangular flow helps to disentangle the ’trivial’ initial 
state fluctuations from the interesting fluctuations that arise due to the critical dynamics. Last but 
not least directed flow is first of all a measure of the initial angular momentum in the heavy ion col¬ 
lision. The beam energy dependence was predicted to show a dip structure when a first order phase 
transition occurs during the evolution. In the following Sections, the beam energy dependence of 
each of these observables will be discussed in more detail. 

3.1 Elliptic flow - a measure of the perfect fluid 

The second Fourier coefficient, the so called elliptic flow, is the response of the system to the 
initial almond shaped overlap region in non-central collisions. Due to fluctuations V 2 is also non¬ 
zero in central collisions, but it is mainly generated due to the overall shape deformation of the 
fireball. Only if the mean free path is small enough the initial spatial anisotropy is converted to the 
corresponding momentum space anisotropy. The beam energy dependence of elliptic flow is quite 
interesting. At very low energies, the spectator nuclei are blocking the interaction region and the so 
called ’squeeze-out’ leads to negative elliptic flow values with respect to the reaction plane. In this 
region the nuclear interactions are important and a hadron transport approach including mean fields 
can describe the elliptic flow rather well (see Fig. 1, left). At high beam energies the elliptic flow 
turns positive and grows as a function of beam energy. In the intermediate region around 10-40A 
GeV the hadron transport approach reaches the right "ball park" values while the underestimation 
at high beam energies is a sign for the importance of partonic interactions starting at around 160A 
GeV. 
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Fig. 1 (right) shows the response function V 2 /e 2 as a function of the charged particle density. 
The green lines indicate the expectations from ideal hydrodynamic calculations and that this line 
meets the experimental data at = 200 GeV was the basis for the claim that a perfect fluid has 
been created. The black lines show various UrQMD hybrid calculations to demonstrate that the 
generic behaviour of this curve can be reproduced once the non-equilibrium stages of the reaction 
are treated properly. 

To investigate how this energy dependence of elliptic flow is generated in the hybrid approach 
the contribution to the final integrated elliptic flow of charged particles has been calculated and the 
percentage is shown in Fig. 2 (left). At lower energies, the hydrodynamic evolution lasts only a 
few fm/c, therefore almost all the elliptic flow is generated by the hadronic transport approach. At 
higher energies the hydrodynamic evolution is necessary to reach the high elliptic flow values as 
discussed above. Fig. 2 (right) shows that the contribution from the late hadronic rescattering is 
constant at about 10 % throughout the whole energy regime investigated. 

Fig. 3 (left) shows the transverse momentum dependent elliptic flow of pions. NA49 mea¬ 
surements are compared to hybrid and pure hadronic transport calculations and again it is clearly 
visible that at 40 AGeV the hydrodynamic evolution does not play a role while at 160 AGeV the 
hydrodynamic evolution is crucial to reproduce the elliptic flow. The dependence on initial state 



Pt IGeV] 



b[fm] 


Figure 3: (Color online) Left: Elliptic flow of pions in mid-central (b=5-9 fm) Pb-rPb collisions at Ei^h = 
4QA GeV and Eiab = 160A GeV. The full and dashed black lines depict the hybrid model calculation, while 
the pure transport calculation is shown as the black dotted line. The colored symbols display experimental 
data obtained with different measurement methods by NA49 [14]. Fig. taken from [12]. Right: Centrality 
dependence of elliptic flow of charged particles at midrapidity (|y| < 0.5) for Pb-rPb collisions at Eiab = 40A 
GeV (a) and Eiab = 160A GeV (b). The horizontal lines indicate the results for averaged initial conditions 
using the hadron gas EoS (blue full line) and the bag model EoS (black dotted line) while the symbols (full 
circles for HG-EoS and open squares for BM-EoS) depict the results for the event-by-event calculation. Fig. 
taken from [15]. 
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Figure 4: (Color online) Left: V 3 {EP} compared with initial state fluctuations’ contribution to V 2 , Oia = 

5 (v 2 {EP}^ — V 2 {RP}^) (squares). Right: Scaled flow coefficients V 2 /{£ 2 ) and vjKej,) with respect to the 
average total hydro duration for impact parameter ranges b — Q — 3A fm, 6.7 — 8.2 fm and 8.2 — 9.4 fm. Figs 
taken from [13]. 

fluctuations and the equation of state has been studied in Fig. 3 (right), where the centrality de¬ 
pendence of elliptic flow of charged particles is shown. There is no difference visible between 
the fluctuating events (symbols) and the averaged smooth initial conditions (horizontal lines) as 
expected for elliptic flow. The difference between the bag model equation of state with a strong 
first order phase transition and the hadron gas equation of state is not observable either. 

3.2 Triangular flow - a measure of fluctuations 

Since 2010 event by event fluctuations in the initial state have been studied as the source of 
higher odd anisotropic flow coefficients [16, 17]. These odd coefficients average to zero without 
event by event fluctuations. To demonstrate that V 3 is really sensitive to fluctuations the beam 
energy dependence of triangular flow of charged particles is compared to the dynamical fluctuations 
of elliptic flow avj in Fig. 4 (left). In central collisions both have roughly the same value whereas in 
non-central collisions they have a similar magnitude only at high energies. At lower beam energies 
triangular flow decreases. 

If non-zero triangular flow is measured at lower beam energies and it is generated by initial 
spatial fluctuations it serves as a nice probe to disentangle these ’trivial’ fluctuations from more 
interesting fluctuations due to the phase transition or critical endpoint. This option relies of course 
on the fact that a finite triangular flow needs to be present. As Fig. 4 (right) shows, V 3 scales for 
different energies and centralities with the duration of the hydrodynamic evolution in the hybrid 
approach. The circles show that for elliptic flow the influence of the hadronic non-equilibrium 
evolution is much larger. Therefore, triangular flow is much more sensitive to a finite viscosity and 
might even vanish at lower beam energies. On the positive side, this disappearance of triangular 
flow can be seen as a more direct evidence that the perfect fluid QGP does not persist for a long 
enough time. 

3.3 Directed flow - a measure of the phase transition 

Last but not least, let us come to one more anisotropic flow observable. The directed flow 
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Figure 5: [Color online] Left: Slope of v\ of protons and pions around mid rapidity extracted from the ideal 
1-fluid calculations with a bag model and crossover EoS. For particle production we applied a Cooper-Frye 
prescription on a iso-chronous hypersurface. Right: Slope of v\ of protons and pions around mid rapidity 
extracted from the ideal 1-fluid calculations with a bag model and crossover EoS. For particle production we 
applied a Cooper-Frye prescription on a iso-energy density hypersurface. Figs, taken from [18]. 


measures the stream of particles within the reaction plane and has opposite sign in the forward and 
backward hemisphere. Nowadays, a rapidity-even vi dipole component has also been observed, but 
here we concentrate on the traditional rapidity odd measure. To plot the beam energy dependence 
one extracts the slope of directed flow at midrapidity. 

Based on the AGS and NA49 measurements there have been predicitions within hydrody¬ 
namic calculations that the slope of the directed flow of protons has a dip as a function of beam 
energy, signaling the phase transition to the quark gluon plasma [19, 20]. In light of the new STAR 
measurements from the RHIC beam energy scan program [21], this old observable has been recal¬ 
culated with modern theoretical techniques. In Fig. 5 a one fluid calculation has been performed. 
This resembles the previous predictions enhanced by actual particle sampling and a v\ calculation 
that can be compared to experimental measurements. On the left hand side, the rather unphysical 
scenario of a constant time freeze-out is shown, whereas the right hand side corresponds to a more 
realistic iso-energy density transition scenario. The two figures demonstrate, that the freeze-out 
dynamics has a significant effect on the structures in the excitation function of the slope of the 
directed flow of protons and pions. The large difference between the two different scenarios - with 
and without first order phase transition - disappears. 

Figs. 6 and 7 both show on the left hand side the directed flow as a function of rapidity at 
two different beam energies and on the right hand side the slope around mid rapidity that has been 
extracted to display the full beam energy dependence. Especially for the proton directed flow it is 
clear that the signal is very small around zero. Even though all the different calculations capture 
the positive slope it is obvious that pure one fluid calculations lead to too much directed flow, while 
surprisingly the pure hadron transport shows the best agreement with experimental data at this 
point. 

Eor pions the negative slope arises because of a shadowing effect. The newly produced parti¬ 
cles are blocked by the spectators and therefore the slope of directed flow is negative at low beam 
energies. Again the fluid calculation overestimates the flow at low beam energies and for pions the 
hybrid UrQMD approach seems to describe the measurements rather well. More important than 
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Figure 6: [Color online] Left: Comparison of proton vi (y) for the different model applied. Right: Slope of 
vi of protons and anti-protons around mid rapidity extracted from the hybrid model calculations with a bag 
model and crossover EoS. We compare with standard UrQMD and experimental data.Figs. taken from [18]. 



Figure 7: [Color online] Left: Comparison of pion vi(y) for the different model applied. Right: Slope of 
vi of negatively charged pions around mid rapidity extracted from the hybrid model calculations with a bag 
model and crossover EoS. We compare with standard UrQMD and experimental data.Eigs. taken from [18]. 


a quantitative comparison to experimental data is the conclusion that the hybrid calculations with 
different assumptions about the order of the phase transition do not show any sensitivity on these 
observables. Detailed measurements of the centrality dependence of directed flow are going to al¬ 
low for better conclusions on the stopping mechanism, the nuclear mean fields and the equation of 
state that are all crucial to understand the dynamics of heavy ion collisions at low beam energies. 

4. Summary and Conclusions 

Collective flow is one of the main observables to understand the properties of the quark gluon 
plasma and the dynamics of heavy ion reactions. Especially the beam energy dependence of 
anisotropic flow provides high potential for insights on the equation of state and the dependence of 
transport coefficients on the baryo-chemical potential. To connect the final state observables to the 
quantities of interest, hybrid approaches based on relativistic hydrodynamics and non-equilibrium 
transport provide a realistic description of the dynamics. The main conclusions are that the small 
change of elliptic flow as a function of the beam energy can be explained by the transport dynamics 
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that gains importance at lower beam energies and replaces the diminished hydrodynamic evolution. 
Triangular flow is much more sensitive to the viscosity and cannot be generated by hadronic trans¬ 
port alone. The beam energy dependence of the slope of the directed flow around midrapidity has 
been studied and even though the old predictions are reproduced, there is no difference between 
a modern hybrid calculation with and without a strong first order phase transition. To reach final 
conclusions on this subject, more detailed studies including nuclear potentials and an improved 
treatment of the baryon stopping are necessary. 
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